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1. INTRODUCTION 

High performance motor drives are necessary parts of industrial applications. It has special 
characteristics such as good dynamic speed command tracking and load regulating response [1]-[4]. The 
methods of control are much simpler and less expensive than those of alternating current motors [5]. DC 
motors are widely used in telecommunication equipment power supplies, stair lifts, household appliances, 
certain low-cost automotive applications and as backup systems in industrial applications [6]-[8]. The 
present 4 quadrant control strategy for DC motors uses the class E chopper which has four active switches 
and requires a dead time between the two switching states [9], [10]. This paper proposes to replace the class 
E chopper with a semi-Z source inverter topology. The proposed topology forms a topological circuit 
consisting of three inductors, three capacitors, a resistor and two switches apart from power supply and back 
emf [11]—{13]. 

DC motor speed is directly proportional to armature voltage, by adjusting the armature voltage 
speed can be varied desirably [14]. Also, the DC motor is capable of running in both directions, which means 
four quadrant operation is possible [15], [16]. So, it is important to design a proper controller to control the 
speed of the DC motor [17]. However, due to lack of an appropriate mathematical model, a proper system 
using the alternative approach could not be designed with as much success as the original 
approach [18]-[20]. Nonetheless, a proper closed loop system using the originally proposed approach was 
successfully designed. 

In order to bridge the gap in area of DC motor drive system this paper attempts the: i) dynamics of 
semi-Z source converter [21], ii) four quadrant operations of semi-Z-source inverter (ZSI) with resistor, 
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inductor, battery load [22], iii) analysis of semi-Z source inverter, iv) designing of closed loop speed control 
of the system. 


2. PROPOSED SEMI-Z SOURCE CONVERTER 

The circuit shown in Figure 1 is the circuit taken up from the base paper and slightly modified. The 
original circuit uses a 2-pole Inductor and Capacitor filter on the output end. This is because the semi-Z 
source inverter was originally developed to be used for single phase Photovoltaic applications which require 
the output voltage ripple to be maintained within certain limits. This is dependent on the nature of the load. 


Figure 1. Semi-Z source inverter circuit 


For a DC machine load, however, the output voltage ripple is not as critical. The machine has a back 
emf which ensures that the armature current ripple is maintained within acceptable limits even with a 
significantly high armature voltage ripple. Hence, the output side LC filter can be removed. The original 
circuit, thus, consists of 2 inductors and 2 capacitors and is a fourth order system. Omitting the LC filter 
reduces the order of the circuit from 4 to 2. This makes the system easier to analyze and work with. 


3. ANALYSIS AND DESIGN 
3.1. Notations and assumptions 

For any variable, the large signal value is the sum of its steady state and small signal value: i) lower 
case denotations represent large signal variables. For example, d represents the large signal value of the duty 
cycle [23]-[25], ii) upper case denotations represent the steady state component of that variable. For 
example, D represents the steady state component of the duty cycle, iii) lower case denotations with a hat on 
top represent the small signal component of that variable. For example, d represents the small signal 
component of the duty cycle, and iv) components with a dot on top represent time derivatives. For example, 
d represents the derivative of the large signal value of the duty cycle with respect to time, v) <d> represents 
the average value of d over a switching period Ts. Assumptions made: i) derivatives of steady state value 
components is zero i.e, D = 0, and ii) product of 2 or more small signal components is negligible i.e, dx 
D5=0. 


3.2. Large signal model 

The large signal analysis for mode 1, when switch Sj is on as in Figure 2 (a): 
for inductor, v,= Lia,, for capacitor, i-= Cvç, and, Vo= Vyn. These equations can be arranged in state space 
form as shown. 


x= A,x + Bu 
y=Cx + Eu 


Similarly, we get the following equations for switching mode 2, when S2 is on as in Figure 2 (b), for 
inductor, v,= LiA= Vm, for capacitor, ic= Cvc= 7 and v= Vc. 


Arranging in state space form, 


x = Áx + Bou 
y = Cx + Eu 


Average values are given as: 
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A=A,d + A,(1—d);B = Bid + B,(1-— d); 
C =C,d + C,(1—d);E = E,d + E,(1-d) 


Thus, we get the final large signal state space model. 


x = Ax + Byu 
y=Cx + Eu (1) 


(b) 


Figure 2. Switching modes in (a) Sı ON and (b) S2 ON 


3.1. Steady state analysis 
The state space average model obtained is not linear. To get a linear system, we need to linearize the 
system around a steady state condition. 


Thus, x=0=A,x + Bu (2) 
The steady state current and voltage of L and C from (2) are given as: 


Ve = —Vim(1-D) & l, = Kea- 


from (1) we see that duty ratio [4] is, Me at 
In 


3.2. Small signal analysis 
Small signal perturbations are added to the large signal model, that is +2 , upon simplification we 
get, 


Ê = AR + Bû + [(4; — A,)X + (B; — By)U]d 

The last 2 elements of the right-hand side can be combined to form a single element such that, 
Ê = AR +[B; (Ay — A2)X + (B1 — BU] inew 

Here, 

If [B ; (Ay — A2)X + (By — B2)U = Brew 


Then, £ = AX + Brewthew 
Similarly, adding perturbations in the output equation, that is Y + ĵ, upon simplification we get, 


J = C + EA + [(C, — C2)X + (E, — E,)U]d 
Finally, Ê = Cê + Enewtnew 
where, [E ; (C, — C2)X + (E, — E3)U = Eney 
Thus, we get the small signal model, 


È = A2 + Bnew bnew (3) 


and, 
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y = CX + Enewtinew (4) 


Applying Laplace transform to the small signal input in (3) and output in (4), and then substituting for 
X(s)we get, 


Y(s) =[C(sI — A) Bnew + Enew]&new(S) 


4. SIMULATION STUDY 

The Semi-Z source inverter is designed considering the output which is intended. The main two 
components are the inductor and capacitor. Using the current and voltage ripple equations, the design 
equations of L and C are as follows: 


Vin (1-D) Ts Vin (1-D)2 Ts 
=———— & C= ~ 
2 Ai, 2 Avc DR 


L 


Based on these equations, the value of inductance is calculated for 4iz-2 A (allowing 2 amperes of 
ripple in inductor current) and the value of capacitance is calculated for Avc-0.125 V (allowing 0.125 volt of 
ripple in capacitor voltage), both at 0.38 duty cycle for continuous conduction. Duty cycle has been limited in 
range 0.38 to 0.75 for proper switching due to asymmetry at extreme points. 

Based on the simulation profile data mentioned in the small signal transfer function of semi-Z 
source inverter is calculated at duty cycle D=0.5, and is given as: 


P(S) _ —1.438s + 7188 
d(s)  0.3125s? + 3.125s + 3.125e04 


From the transfer function it can be clearly seen that this semi-Z source inverter topology has a right half zero 
which results in a non-minimum phase system. 

The speed is taken as feedback and is compared with reference speed, the generated error signal 
passes through the outer loop Proportional Integral controller which generates current proportional to 
reference voltage which is then compared further with feedback voltage in the inner loop. This generated 
error signal passes through the inner loop proportional integral derivative (PID) controller and generates a 
current proportional to duty ratio. Saturation block is used in the outer loop to limit voltage in range -150 to 
150 and saturation block in inner loop is used to limit duty ratio in range 0.38 to 0.75. Pulse width 
modulation (PWM) pattern is generated for the generated error signal by PWM block which is given to 
complementary switches S; and S2. In this way this particular closed loop converter works. 

According to the study by the authors, it is a non-minimum phase system having a right half zero, so 
it exhibits an undershoot at start. Due to this undershoot, error increases and while trying to correct this error, 
the controller increases duty cycle due to which it reaches peak negative value and then at this point it finally 
corrects itself, which is not acceptable. To overcome this, integral gain needs to be increased, which makes 
the response sluggish, which is acceptable over previous cases, which is shown in Figure 3. 


Speed & Arm. current (A) 


o 0.5 1 1.5 2 253 3 3.5 4 4.5 5 
Time (S) 


Figure 3. Closed loop systems with small integral controller exhibiting undershoot at start 
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4.1. Tuning of PID controller 

The classical controller design methods are used because physical systems are easier to implement 
and design as compared to one designed using modern control theory. These controllers are preferred in most 
industrial applications. For inner loop Ziegler Nichols sustained oscillations method is used, whereas for 
outer loop Chien, Hrones and Reswick Open loop tuning technique is used. 

The controller transfer function is given as: 


1 
G.(s)=Kp(1 +t Pes) 
Both the methods are discussed. 


4.1.1. Ziegler Nichols closed loop tuning or sustained oscillations method 

This method of tuning is limited to tuning processes that cannot be applied on open loop systems. It 
allows us to find critical gain Kc; and critical period of oscillations Pc. A closed loop system with unity 
feedback is used in this method [15], [18]. The integral and derivative terms are made zero, with only 
proportional gain (K). Initially K is set at low value and is gradually increased till the verge of instability, 
when the system produces sustained oscillations, it is shown in Figure 4. Gain at which the system produces 
sustained oscillation is called critical gain Kcr and the period of oscillation is Pcr. Based on both the 
parameters PID controller gains are calculated by Table 1. Based on the performance of the system, values of 
critical gain and critical period are found as Kc=2.0201 and Pc.=0.016 s respectively, which is shown in 
Table 2. Sustained oscillation with period Pc; for gain Kc;is shown in Figure 5. 

So based on the calculated values of Kp, Ti, Ta controller transfer function is given as: 


151.512 


G(s) = (1.21244 4240-35) 


Output 


Speed Voltage 
[TH | —1.438s + 7188 = 


0.31255 + 3.125s + 3.125e04 


Converter Tranfer Function 


Figure 4. Closed loop inner voltage loop with proportional gain 


Table 1. Ziegler Nichols closed loop calculation of Table 2. Calculated values of Kp , Ti & Ta for PID 


Kp, Ti & Ta controller 
Type of Controller Kp Ti Ta Type of Controller Kp Ti Ta 
PID 0.6Kce, 0.5Pc,  0.125Pc: PID 1.212 0.008 0.002 
35 ~~ 
30 
25 
= 
Z 20 
wo 
Fi 15 
$ 10 
5 
o 
0 0.5 1 1.5 2 2.5 3 3.5 4 45 5 
Time (S) 


Figure 5. Sustained oscillation with period Pc; for gain Kcr 
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4.1.2. Chien, Hrones and Reswick (CHR) open loop tuning technique 

This is a modified version of Ziegler Nichols open loop process reaction method which provides a 
better way of selection of controllers for control applications. The process of finding controller gains for 
setpoint regulation with 0% overshoot is summarized is shown in Figure 6. 


85.3 % 


Amplitude (A) 


35.3 % 


t1 t2 Time (S) 


Figure 6. Step response of open loop process for CHR technique 


A=1.3 t1— 0.2 t2 
B=0.67 (ti — t2) 
K=Output at steady state/input at steady state 
Through the given standard formula, the PID controller gains can be calculated is shown in Table 3. 
Based on the results tı and tz are found as tj=1.45 s and t2=16 s respectively, with these values, parameters for 
CHR method is shown in Table 4. And the Figure 7 shows Open loop simulation of outer speed loop for 
CHR technique. Step response of outer open loop is shown in Figure 8. 
A=0.09 
B=8.375 
K=4.1 
So based on the calculated values of Kp, T;, Tacontroller transfer function is given as: 


1.626 


G.(s) = (13.618+——+ 51.025s) 


Table 3. CHR open loop calculation of Kp, Ti , Ta Table 4. Calculated values of Kp, Ti & Ta 
Controller Kp Ti Ta Controller Kp T; Ta 
PID 0.6/K _ 2 05 PID 13.618 8.375 4.188 


Speed 
ig -1.4385 +7188 
E 0.31252 +3.125s + 3.125204 2+ 50.1255 + 6.25 


Converter Tranfer Function DC Motor Transfer Function 


Inner Voltage Loop feedback 


Figure 7. Open loop simulation of outer speed loop for CHR technique 


Voltage (V) 


0 01 02 03 0.4 05 06 07 08 0.9 
Time (S) 


Figure 8. Step response of outer open loop 
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4.2. Tuning of PID controller - manual fine tuning 

The PID controllers tuned using classical control design methods involve design approximations, so 
the fine experimental on-site tuning is necessary to yield accurate results. The following are the controller 
gain values as obtained by manual tuning and the dynamics of change in controller gains, it is shown in 
Table 5. 

The tuning of the inner loop is done first followed by the outer loop. The derivative controller is not 
used in the control system unless necessary, that is why most of the system employs PI controllers only. This 
is because when a sudden input is applied to the system, it compensates the output fast which in turn affects 
the controller to allow huge steady state errors in the long term. But in this system there is an undershoot in 
the negative direction of control, so derivative action needs to be taken into account to improve the transient 
response of the system. It should be kept in mind that increasing the value of Kp beyond range slows down 
the response. 

Therefore, based on all the analysis done so far, the final controller transfer functions are as follows: 


PID Inner loop: Ge(s) = 1.2 + = + (2.424e — 3s) 


PID Outer loop: G(s) = 20 + = + (10s) 


Table 5. Controller gains and system dynamics 


Gain Value Increase Decrease 
Inner voltage loop 
Kp 1.2 Little overshoot in current No significant chance (slow response) 
Kı 160 Overshoot Undershoot at start 
Kp 2.4e-3 Slows down response Undershoot at start 
Outer speed loop 
Kp 20 Overshoot in current Sustained oscillations 
Kı 1.6 Overshoot in speed Ripples in current 
Kp 10 Slow response Undershoot for negative speeds 


5. RESULTS AND DISCUSSIONS 
5.1. Open loop simulation results 

The following figures show the open loop simulation output done using MATLAB/Simulink for all 
the four quadrants with varying duty ratios. Initially the duty cycle is set to 70% at t=0 s, then 50% at t=2.6 s, 
then 40% at t=5 s and finally 60% at t=7.6 s. For this change the duty cycle speed changes from 0 to 62 rad/s 
(forward motoring), then from 62 to -10 rad/s (forward braking), then from -10 to -80 rad/s (reverse 
motoring) and finally from -80 to 32 rad/s (reverse braking) respectively, current is shown in Figure 9 and 
open loop-output voltage and it’s mean is shown in Figure open loop-speed and armature 10. 


Speed & Arm. current (A) 


2.5 3 3.5 4 45 5 
Time (S) 


0 0.5 1 15 


N 


Figure 9. Open loop-speed and armature current 
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Speed & Arm. current (A) 


o 0.5 1 15 2 25 3 3.5 4 45 5 
Time (S) 


Figure 10. Open loop-output voltage and it’s mean 


5.2. Closed loop simulation results 
5.2.1. For varying speed reference 

The simulation results of the closed loop system for varying speed reference with PID controller is 
shown below. The speed reference is set to 50 rad/s at t=0 s, 30 rad/s at t=2.6 s, -50 rad/s at t=5.1 s and 
finally to -30 rad/s at t=7.5 s. Figure 11 shows closed loop-speed and armature current and closed loop-output 
voltage and it’s mean is shown in Figure 12. 


Speed & Arm. current (A) 


o 0.5 1 1.5 2 25 3 3.5 4 45 5 
Time (S) 


Figure 11. Closed loop-speed and armature current 


Speed & Arm. current (A) 


o 0.5 1 15 2 2.5 3 + 4 45 5 
Time (S) 


Figure 12. Closed loop-output voltage and it’s mean 
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5.2.2. For varying load torques 

Shown below are the simulation results of the closed loop system for varying Load Torque with PID 
controller. The load Torque is varied from 0 to 10 N.m at t=2.5 s, from 10 to 50 N.m at t=5 s and finally from 
50 to 30 N.m at t=7.5 s. It can be seen that speed remains constant for different load torques. Closed loop- 
speed and armature current is shown in Figure 13 and closed loop-load and electromagnetic torque is shown 
in Figure 14. 


Speed & Arm. current (A) 


2.5 3 35 4 4.5 5 
Time (S) 


o 0.5 1 LS 


N 


Figure 13. Closed loop-speed and armature current 


Speed & Arm. current (A) 


o 0.5 1 15 2 25 3 3.5 4 45 5 
Time (Ss) 


Figure 14. Closed loop-load and electromagnetic torque 


6. CONCLUSION 

In this paper a novel approach for the 4-quadrant speed control of a DC motor using a semi-Z source 
inverter was introduced. The steady state and small signal analysis of the semi-Z source inverter topology 
were discussed. The topology was first tested for open loop operation and it was found that it meets the 
requirements for the intended operation. The results of the analysis of the topology were used to design a 
closed loop system for output voltage control in both polarities. A DC machine load was connected to the 
output terminals and another closed loop control system was designed for speed control of the machine in all 
four quadrants. The designed system was simulated using MATLAB/Simulink and the results were discussed. 
The designed system successfully regulated the speed of the DC machine in all four quadrants and managed 
to maintain a constant speed even in varying torque conditions. Thus, it was found that the overall system 
works as intended. The designed system successfully meets all the performance requirements as well as the 
proposed objectives. 
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